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ABSTRACT 



Q . We present infrared (IR) spectrophotometry {R ~ 180) of three Herbig Ae/Be 

stars surrounded by possible protoplanetary disks: HD 150193, HD100546 and 
HD 179218. We construct a mid-IR spectral energy distributions (SED) for each 
object by using 7.6 — 13.2 /im HIFOGS spectra, 2.4 - 45 /im spectrophotometry 



c/3 . from the ISO SWS, the 12, 25, 60, and 100 /xm photometric points from IRAS, 

^ ' and for HD 179218, photometric bolometric data points from the Mt. Lemmon 

Observing Facility. The SEDs are modeled by using an expanded version of the 
^ ! Chiang & Goldreich (1997) two-layer, radiative and hydrostatic equilibrium, pas- 

sive disk. This expanded version includes the emission from Mg-pure crystalline 
olivine (forsterite) grains in the disk surface layer. Each of the three objects 
studied vary in the amount of crystals evident from their spectrophotometry. 
HD 150193 contains no crystals while HD 100546 and HD 179218 respectively 
show evidence of having crystalline silicates in the surface layers of their disks. 
We find that the inner region of HD100546 has a 37% higher crystalline-to- 
amorphous silicate ratio in its inner disk region (< 5 AU) compared to the outer 
disk region, while the inner disk region of HD 179218 has a 84% higher crystalline- 
to- amorphous silicate ratio in its inner disk region (< 5 AU) compared to the 
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outer region. All three objects are best-fit using a grain size distribution power 

law which falls as ar^-^ . HD 150193 is best-fit by a small disk (~ 5 AU in radius) 
while HD 100546 and HD 179218 arc best fit by larger disks (~ 150 AU in radius). 
Furthermore, HD 100546's disk flares larger than compared to HD 150193 (25% 
more at 5 AU) and HD 179218 (80% more at 5 AU). We discuss the implications 
of our results and compare them with other modeling efforts. 

Subject headings: stars: general 

1. INTRODUCTION 

The ubiquity of accretion disks around pre-main sequence and young main sequence 
stars having the potential to form planetary systems is now well established (Koerner 2001). 
However, unknown is an accurate estimate of the fraction of single stars with disks that have 
produced planetary systems. If not all such systems have planets, why then do only some 
accretion disks form planets and others do not? Theoretical models of particle aggregation 
show that if particles can grow from submicron to mm to cm in size, then the formation of 
planetesimals is possible in the time before the disk dissipates (Habing ct al. 1999). The 
problem remains to understand how grains condense from nebular gasscs, and how relic 
interstellar grains survive and arc modified by their transport in the disk. These pristine 
submicron size grains, be they nebular condensates or relic interstellar grains, then grow into 
mm size particles. However, investigating how grains grow is complicated by the fact that 
most of the grain processing probably occurs in the hot mid-planes of these disks (Bell et al. 
2000), hidden from view by the optically thick disk photospheres. If grains are lofted above 
the photosphere by disk processes such as winds, turbulent convection, or changes in vertical 
structure, the evolution of dust can be investigated by observing the properties of the small 
(< 1 yum) grains in the disk surface layer or atmosphere. 

In particular, observing and modeling the 10 /xm silicate resonances in nearby (d < 
500 pc), young (age < 5 Myr) Herbig Ae/Be (HAEBE) stars, the sites of possible on-going 
planetary formation (Waters & Waelkens 1998), establishes the silicate mineralogy, Mg con- 
tent of the grains, the degree of crystallinity vs. amorphousness, grain porosity, and grain size 
distributions. These grain properties can be linked to conditions in protoplanetary disks and 
constrain theoretical models of protoplanetary disk evolution. For example, amorphous sili- 
cates commonly detected in solar system comets are considered to be interstellar rehc grains 
(Li & Greenberg 1997) and probably are represented by Glasses Embedded with Metals 
and Sulfides (GEMS) in interplanetary dust particles, IDPs (Brownlee et al. 2000). Con- 
versely, Mg-rich crystaUine silicates are thought to be either pristine solar nebula condensates 
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that condensed at ~ 1450 K or amorphous silicates annealed into crystals at temperatures 
> 1000 K in the inner hot zones of the accretion disk. If these silicate crystals are formed 
in the hot inner zones of the disk, then the crystals must migrate to large radial distances 
in order to be incorporated into comets and outer nebula protoplanetesimals. However, 
the hot temperatures required to form the silicate crystals could have occurred in the early 
phases (~ 300, 000 yrs) of the solar nebula in the disk (Brocklee-Morvan et al. 2002) or 
in nebula shocks in the 5 — 10 AU region, obviating the need for radial transport of the 
crystals to comet forming zones. Indeed, Mg-rich sihcate crystals within IDPs are detected 
through cometary IR spectra (Wooden et al. 2000a; Wooden, Marker, & Woodward 2000b). 
Crystalline silicates also have been spectroscopically identified in pre-main sequence HAEBE 
stars including HD 100546 (Malfait et al. 1998) and in P Pic (Pantin, Lagage, & Artymowicz 
1997), a debris disk system continually repopulatcd with dust grains via comet-comet colli- 
sions and cometary accretion events. The HAEBE stars (the high mass [2 — 10 Mq] analogs 
of T Tauri stars [~ 1 M0]) are the probable progenitors of the /3-Pic stars since both have 
a common occurrence of 'exocometary activity,' are similar in spectral type, and have an 
apparent overlap in Hipparcos ages (Crifo et al. 1997). 

In this work, we examine the thermal emission from three Herbig Ae/Be stars of similar 
spectral type: HD 150193 (AlV), HD 100546 (B9V), and HD 179218 (B9). The stellar ages 
of the objects are: 2 x 10^ yrs, 10 x 10^ yrs, and 0.5 x 10^ yrs for HD 150193, HD 100546 
and HD 179218, respectively (van den Ancker et al. 1998). We chose these objects for study 
because they are isolated (not belonging to any know star forming region), have little to no 
active accretion, possess possible processing circumstellar disks, and have varying degrees of 
silicate crystallinity (Meeus et al. 2001). We assemble and model SEDs for each system using 
a passive reprocessing circumstellar disk to constrain the amount of crystalline silicates in 
the inner and outer parts of the disk, and flaring of the disk in each of the systems. 



2. OBSERVATIONS 

Spectrophotometry of HD 150193, HD 100546, and HD 179218 was obtained using the 

NASA Ames Research Center High efficiency Faint Object Grating Spectrometer (HI- 
FOGS: Witteborn et al. (1991)). The spectrum of HD 179218 was obtained on the 2.34-m 
telescope at the Wyoming Infrared Observatory (WIRO) using 32-bit FORTH telescope 
software. The spectra of HD 150193 and HD 100546 were obtained on the 4-m Blanco 
telescope at the National Optical Astronomical Observatories Cerro Tololo Inter- American 
Observatory (CTIO) in Chile. Table 1 summarizes the HIFOGS observations used in this 
work. 



-4- 



The HIFOGS is a grating spectrometer with a resolution of AA ~ 0.045 //m per detector, 
yielding a resolving power oi R = A/AA ~ 350 — 180 covering the 7.5 — 13.4 /xm spectral 
region. A 3" diameter circular entrance aperture was used for each observation at each 
observing facility. This aperture corresponds to a radial size of 448 AU, 309 AU, and 732 AU 
for HD 150193, HD 100546 and HD 179218, respectively. Standard IR observing techniques 
of chopping and nodding were employed, with a chop frequency of 7 Hz and with a chop 
throw of 30" in the El-W direction. Absolute cahbrated flux spectra of the standard stars 
used in the data reduction were taken from the UKIRT CGS 3 and HIFOGS measurements 
compiled by Cohen et al. (1996). Ratios of two or more measurements for each standard 
at different air masses were fitted to ratios of atmospheric transmissions (using ATRAN 
software. Lord (1993)), for a precipitable water vapor of 4.5 mm for both WIRO and CTIO to 
verify the wavelength calibration and the atmospheric transmission corrections as functions 
of airmass. Each HIFOGS spectrum was flux calibrated by dividing the measured spectrum 
by the standard star, multiplying by the corresponding star's flux spectrum (Cohen et al. 
1996), and multiplying by the inverse ratio of computed atmospheric transmission spectra. 
All observations presented here were obtained under photometric conditions. The spectra 
are shown in Fig. 1. 

Complete Nyquist sampling of the 10 /xm spectral signatures requires two HIFOGS grat- 
ing settings. We were able to achieve Nyquist sampling for the object HD 179218. The other 
two objects, HD 150193 and HD 100546, were observed at a single grating setting. Obser- 
vations of HD 179218 on 1998 June 23 UT, June 27 UT, and June 30 UT were obtained at 
the first grating setting that spanned 7.811 — 13.356 iim. The 1998 June 24 UT and June 
25 UT observations were obtained at a second grating setting that was shifted such that the 
wavelength interval of the detectors was half way between the wavelength interval of the de- 
tectors in the flrst grating setting. The second grating setting spanned a wavelength range of 
7.519 — 13.122 /xm. The spectra from each observational epoch are shown in Fig. 1. Spectra 
for each grating setting were scaled to the spectrum obtained during the best photometric 
conditions (1998 June 23 UT and 1998 June 25 UT for the first and second grating settings, 
respectively) and combined by computing a statistically weighted average. The scaling fac- 
tors for each spectrum are hsted in Table 1. The resultant HIFOGS spectra for HD 179218 
at each grating setting were subsequently interlaced to produce the final spectrum (spanning 
7.528 — 13.356 ^m). Therefore, for HD 179218, there are approximately twice as many data 
points between 7.5 fim and 13.4 /im than for HD 100546 and HD 150193. No scaling factor 
was used to interlace the spectra from the two grating settings. 

We also obtained 1.2 - 23.0 /xm IR broad- and narrow-band photometric observations 
of HD 179218 on 2001 April 30.42 UT using the University of Minnesota (UM) Mount 
Lemmon Observing Facility (MLOF) 1.52 m telescope and the UM multi-filter bolometer 
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(Manner et al. 1990). Observations where conducted using 9.33" circular aperture, a 29" 
N-S nod, and a 7.5 Hz chop-frequency. In addition to standard near-IR broad-band filters, 
the bolometer filter set also included the six Infrared Telescope Facility (IRTF) narrow-band 
7-13 fim "silicate" filters (Tokunaga et al. 1986). Multiple observations of the photometric 
standard star a Boo (= HD 124897) at a variety of air masses during the night established 
the photometric zero point over the entire spectral range of the bolometer. The photometric 
systems, magnitude scales, and absolute flux calibrations for the bolometer are given by 
Mason et al. (2001) and Gehrz (1997a,b). Standard extinction values of 0.2 magnitudes per 
air mass (mag AM~^) were used for fllters [J] — [L] and 7.9-12.6 //m, 0.3 mag AM"-*^ for 
fllter [M], and 0.5 mag am~^ at 18.3 fim. The use of standard extinction values introduces 
no more than a few percent uncertainty in the derived photometry, as observations of HD 
179218 where obtained at < 1.2 air masses. Observed magnitudes corrected for atmospheric 
extinction, are summarized in Tables 2. 



2.1. Assembly of Data Sets for SEDs 

Data from the European Space Agency Infrared Space Observatory (ISO) Short- Wavelength 
Spectrometer (SWS) and the Infrared Astronomical Observatory (IRAS) are used to expand 
the wavelength coverage of the spectral energy distribution (SED) of each of our three 
HAEBEs. The ISO SWS observations were conducted using the SWSOl scanning mode (As- 
tronomical Observing Template, referred to as AOTs) covering the entire 2.4-45 /xm SWS 
wavelength range. Each object has been observed at a different scanning speed; in particu- 
lar, data have been taken at "speed 4" (HD100546), "speed 3" (HD179218), and "speed 1" 
(HD150193), corresponding to raw resolution of approximately A/AA ~ 1400, 800, and 400, 
respectively. 

Since the HAEBEs were faint and the standard pipeline process did not adequately 
apply a dark subtraction to the data points, we have reprocessed the data from a raw 
level to a fully processed stage using the Observer's Spectral Interactive Analysis Package 
(OSIA v2.0) and the ISO Spectral Analysis Package (ISAP v2.1). Each of the 12 bands into 
which the SWSOl is divided have been processed separately; dark cm^rcnt subtraction, scan 
direction matching and fiat-fielding have been applied interactively to each band. Using 
ISAP, we sigma clipped outliers and averaged the data of all 12 detectors for each AOT 
band, retaining the instrument resolution. For those bands affected by fringes, we applied 
the defringe routines incorporated in ISAP. 

For extended sources a fiux discontinuity may occur between band limits as a conse- 
quence of a change in aperture size at some band edges. Our three objects show very little 
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flux jumps at the detector band limits, and in most cases, we did not have to normahze 
the segments to their neighbors when we combined the 12 spectral segments into a single 
spectrum. For HD 179218 we use an additive term to normalize the last two segments 
(19.5 — 45 /xm) to the flux in the previous band. Since we have ISO SWS data from objects 
that have been taken at different scanning speeds and, consequently, have different resolu- 
tions, we have produced the final spectrum for each object rebinning the data to a bin size 
that corresponds to a spectral resolution of A/AA ~ 400. 

The ISO SWS spectra are scaled to the HIFOGS spectrophotometry using the following 
procedure: 1) the resolution of the ISO SWS spectrum is degraded to match the wavelength 
scale of the HIFOGS spectrum; 2) a scaling factor is calculated for each data point (ap- 
proximately 30 points for HD 100546 and HD 150193, and 60 points for HD 179218) in the 
high S/N region between 8 and 10 //m; 3) a statistically weighted average of the scahng 
factors is computed; and 4) the averaged scaling factor is applied to the entire original, un- 
degraded ISO SWS spectrum. The scaling factor of the ISO SWS spectra to the HIFOGS 
spectrophotometry for each object is listed in Table 3. 

In addition to spectra from ISO SWS, we use IRAS flux values from Oudmaijer et al. 
(1992) for each of the three objects. We apply the method outhned in the IRAS Explanatory 
Supplement (1988) to produce color corrected, monochromatic flux points. We assume the 
IRAS fluxes are larger due to the larger IRAS aperture. First, the 12 /im IRAS data point 
is scaled to the HIFOGS spectrum. Then, the resulting scaling factor is applied to each of 
the other three IRAS data points (Table 3). 

Finally, we added the IR photometric points from MLOF to the SED of HD 179218. 
The near-IR points were corrected for interstellar reddening by assuming the extinction law 
of Mathis (2000) and the published extinction for HD 179218 of = 1.6. The scaling 
factor applied to the MLOF data to match the HIFOGS spectrum of HD 179218 is listed in 
Table 3. The flnal assembled SEDs for all three objects are shown in Fig. 2. 

The factors needed to scale the ISO SWS data to the HIFOGS data arise due to a 
difference in aperture size and distance to the object. The effective aperture size of ISO SWS 
is 9.44", and for HIFOGS is 3". Based on millimeter studies, the diameter of the disk around 
HD 179218 is 230 - 450 AU (Mannings & Sargent 2000) and for HD 150193 is ~ 250 AU 
(Mannings & Sargent 1997). The distance to HD 179218 (based on Hipparcos data) is 
244 pc and for HD 150193 is 149 pc. This results in projected diameters of the ISO SWS and 
HIFOGS apertures on HD 179218 of 2302 AU and 732 AU. Similarly, the projected diameters 
on HD 150193 are 1409 AU and 448 AU. Therefore, the apertiues contain the millimeter 
measured disk sizes resulting in ISO SWS to HIFOGS scaling factors close to 1 for these 
two objects. Conversely, the distance of HD 100546 is 103 pc yielding a projected aperture 
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diameter from ISO SWS and HIFOGS as 973 AU and 309 AU, respectively. Although there 
are no millimeter measurements of HD100546, HST Chronographic Imaging reveals a disk 
approximately 760 - 1030 AU in diameter (Augereau, Lagrange, Mouillet, & Menard 2001), 
(Grady ct al. 2001). Therefore, significant mid-lR fiux may be contained in the ISO SWS 
aperture which would account for the comparatively large scaling factor needed to scale ISO 
SWS to HIFOGS for HD 100546. 



3. Passive Disk Model 

We use a modified version of the radiative, hydrostatic models of passively irradiated 
circumstellar disks developed Chiang & Goldreich (1997) (hereafter CG97) and updated by 
Chiang ct al. (2001) (hereafter COl). The CG97 disk is composed of two parts: an optically 
thin surface layer and an optically thick interior (Fig. 3). The surface layer is directly 
heated by hght from the central star. Half of the reprocessed emission from the surface 
layer escapes into space while the other half radiates onto the disk interior. The emission 
from the disk surface layer heats the cooler disk interior. Since the disk is in hydrostatic 
equilibrium, the disk flares with increasing distance from the central star. The amount of 
disk flaring is influenced by the temperature of the interior. Flaring of the disk results in 
more starlight hitting the disk at larger distances from the central star as compared to an 
otherwise geometrically flat disk. COl updated the CG97 model by adding a range of grain 
sizes and incorporating a simple mineralogy through the use of laboratory determined optical 
constants, including water ice, amorphous olivine, and metallic Fe. Although the two-layered 
disk approach is simple, when compared to more detailed vertical structure models, the CG97 
model was found to be a robust model for computing SEDs from HAEBEs (DuUemond, van 
Zadelhoff, & Natta 2002). 

In this work, we use the updated model of COl with our own enhanced feature of adding 
Mg-pure crystaUine ohvine grains (forsterite) into the optically thin surface layer. We also 
are able to vary the crystalline-to-amorphous silicate ratio in the inner and outer radial 
regions of the disk. The details of the standard disk model are explained in CG97 and COl; 
however, we will highlight the disk and dust parameters here. 

Disk parameters include the disk surface density at 1 AU (So), the outer disk radius 
(ro), and the height of the disk photosphere in units of the gas scale height (H/h). Dust grain 
parameters include grain mineralogy, slope of the grain size distribution (g), maximum grain 
size {ama:,:}- ^ud subllmatiou temperature {Tsub)- The shape of the SED in the region covered 
by our data sets (2 - 45 /im) are largely unaffected by changes in the input parameters. The 
dust parameter which most affects the shape of the resonance features is the slope of the grain 
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size distribution, q (COl). Changes in disk radius, disk surface density, and maximum grain 
size for the interior and the surface layer mostly affect the shape of the SED at wavelengths 
longer than 50 /xm (COl; Creech-Eakman et al. 2002). The one disk parameter that produces 
the most change at wavelengths less than 50 fim is H/h, which can be used as a measure of 
dust settling in the disk (COl). However, as noted in COl, H/h is not a fixed parameter and 
should be calculated self-consistently at each radius in the disk. We follow the procedure 
of COl by keeping H/h constant throughout the disk and using it as an indicator of dust 
settling in the disk. Another unknown parameter expressed in terms of H/h may actually 
represent the settling of dust in the disk. Finally, the relative mass fraction of the minerals 
in the disk and surface layer will influence the SED in the mid-IR. 

To model the emission from dust particles, we use the method of COl and assume that 
the disk is composed of a size distribution of isolated spherical grains of discrete mineralogy; 
i.e., we ignore gas as a source of opacity. The specific mineralogy is addressed in the next 
section (§3.1). The opacity of the grains is computed by using Mie Theory. For grains coated 
with a water ice mantle, Mie- Guttler Theory is used to compute the opacity. The grain size 
distribution in the disk surface layer and in the disk interior are equal, and unchanged in 
each of our models. The central stars are modeled as simple blackbodies. The radius, 
temperature, and luminosity of the central stars of our three objects are taken from the 
compilation by van den Ancker et al. (1998). We use the luminosity to calculate the radius 
of each of the stars using the relation: L*(ergs cm~^ s^^) = AnaRlT^. 

Two final assumptions are made when using the COl model. The first is that we are 
viewing the disk face-on. The stars discussed here have low A^ {A^ = 0.28 for HD 100546; 
Ay = 1.27 for HD 179218; and Ay = 1.6 for HD 150193), therefore, the disk is not inclined 
at such an angle as to significantly block the starhght. CG99 showed that until the disk is 
inclined at such an angle as to block the starhght {i ~ 45°), the SED is mostly unaffected by 
disk inclination. Therefore, we use the approximation that the disks are face on for all three 
of our objects. The second assumption is that any accretion in these systems is negligible and 
does not significantly contribute to the measured emission from these systems. By making 
such an assumption, we are placing an upper limit on the contribution of flux from the disk 
alone on the measured emission from each of these systems. 

Flux from the central star, disk interior and disk surface are co-added to produce the 
model SED. The relative amount of amorphous to crystalline silicates in the disk surface 
layer is adjusted to produce the best fit within a 2a (95% confidence) level (Press et al. 
1992). 
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3.1. Dust Mineralogy 

The mineralogy used by COl is a good representation of the basic minerals extant in 
young disk systems (Malfait et al. 1998; Bouwman et al. 2000). We maintain the mineralogy 
of COl, using amorphous olivine, metallic Fe, and water ice. However, we have expanded 
the mineralogy by including Mg-pure crystalline olivine grains in the disk surface layer. 
Crystalline olivine has been discovered to be a very important component of the dust in 
circumstellar disks and solar system comets. In an attempt to quantify the amount of 
crystalline olivine grains, we make this important step of including them in the COl passive 
disk model. Table 5 lists the minerals used in the modeling and the reference for the indices 
of refraction. 

It is difficult to model the thermal emission from crystalline silicate grains (Yanamandra- 
Fisher & Hanner 1999). Many authors choose to use a form of a continuous distribution of 
ellipsoids while others use direct comparisons from laboratory transmission experiments. In 
this work, we choose to use the methods outlined by Fabian et al. (2001) who calculate 
the emission from ellipsoids by elongating the crystals along one of the three crystalline 
axes. Based on the location of the resonance peaks in the HIFOGS and ISO SWS spectra 
of HD 150193 and HD 179218, we use an axis ratio of 10:1:1 to compute the optical effi- 
ciencies (Qs) of the crystals. The crystalline olivines are not coated with ice since to do so 
requires mixing theory which eliminates the distinctive crystalline resonances. We compute 
the thermal emission from the crystals in the disk surface layer of each object. We make the 
simplified assumption that the crystals are at all radii in the disk surface layer. This assump- 
tion may not be valid since ice should coat the crystals out at larger disk radii. However, the 
measured crystalline resonances seem to originate from bare crystalline grains. Finally, we 
vary the ratio of amorphous to crystalline silicates to find the best model fit to the SEDs. 

We note that in addition to solid-state emission features from silicates, HD 100546 and 
HD 179218 exhibit emission from polycyclic aromatic hydrocarbons (PAHs) (Meeus et al. 
2001). We do not attempt to model the emission from these molecules here, but will examine 
them in a future paper. 

4. RESULTS 

Our modeling results are displayed in Figs. 4 and 5 and summarized in Table 6. To 
produce the best-fit model to the SEDs, HD 100546 and HD 179218 are both modeled with 
a disk 150 AU in radius, while HD 150193 is modeled with a disk 5 AU in radius. The 
best fit SEDs for HD 100546 uses a scale height of H/h = 4, a scale height of H/h = 3 for 
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HD 150193, and a scale height of H/h — 1 for HD 179218. This means that at a radius of 
5 AU the disk around HD 100546 flares about 80% more than the disk around HD 179218, 
and about 25% more than the disk around HD 150193. The current model fits suggests that 
the measured emission from HD 150193 is dominated by warm dust in a flaring disk close 
to the central star, and that the emission from grains in an extended disk (> 5 AU) are 
relatively small compared to the other two objects. We should note that for HD 150193, a 
larger disk of radius 50 AU with a scale height of H/h — 1 can be used to match the contrast 
of the measured silicate resonance features in HD 150193, but the calculated flux at longer 
wavelengths exceeds the measured IRAS data point at 60 jivci by a factor of 4. Therefore, we 
are able to match the long wavelength photometry and the contrast of the sihcate resonance 
features by reducing the radius of the disk and increasing the scale height. 

The quahtative crystaUine content of these three HAEBEs has previously been reported 
by Meeus et al. (2001) with HD 100546 and HD 179218 showing evidence of crystalline 
silicates, while the emission from HD 150193 is dominated by amorphous silicates. For the 
two objects with crystalline silicates, we find that a better fit is produced to the observed 
SEDs of the objects if we use a model in which a higher ratio of crystalline-to-amorphous 
silicates is located in the inner regions of the disk (< 5 AU) compared to the outer regions 
of the disk (5 - 150 AU). The inner region of HD 100546 has 37% higher crystalline-to- 
amorphous silicate ratio compared to the outer region. This is contrary to the findings of 
(Bouwman, de Koter, Dominik, & Waters 2003) (hereafter, BdKDW) who used a spherical 
shell model to calculate a factor of almost 10 higher fraction of crystaUine sihcates in regions 
greater than 10 AU. The inner region of HD 179218 has 84% higher crystalline-to-amorphous 
silicate ratio compared to the outer region. Our model results are consistent with the recent 
findings of van Boekel et al. (2004). van Boekel et al. (2004) observed three Herbig Ac stars 
(HD 142527, HD 163296, and HD 144432) using the Mid-Infrared Interferometric Instrument 
on the Very Large Telescope Interferometer. In all three of their objects, they found evidence 
for a larger fraction of crystalline silicates in closer to the star (1-2 AU) compared to the 
outer region of the disk (2-20 AU). 

Finally, all three of our objects are best fit using a grain size distribution with a slope 
oi q = 3.5 for both the optically thick inner disk and for the optically thin surface layer. 

5. DISCUSSION 

It is difficult to make any statistically significant conclusions about disk evolution from 
modeling only three objects. However, we can make some interesting observations based 
on our results. Meeus et al. (2001) defined two groups of Herbig Ae/Be stars based on the 
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shape of their SEDs. Group I objects exhibit near-IR and far-infrared (far-IR) emission, with 
and without sihcate emission features (denoted Groups la and lb, respectively). Group II 
objects exhibit near-IR emission and sihcate emission features, but much less far-IR emission 
compared to the objects in Group I. HD 179218 and HD 100546 are considered Group la 
objects and HD 150193 is considered a Group II object under this scenario (Meeus et al. 
2001). 

Meeus et al. (2001) suggest a disk geometry interpretation of their groupings. Group 
I objects have a geometrically thin, optically thin inner region with a flaring outer region. 
Group II objects contain an optically thick inner region which puffs up, shielding the opti- 
cally thin, low mass, outer disk region from stellar flux, thereby preventing it from flaring. 
Qualitatively our modeling results match these groupings. HD 150193 is modeled with a 
relatively small (5 AU), low mass disk, which suppresses the amount of emission from ice 
coated silicate grains contributing to the emission around 45 //m. HD 179218 and HD 100546 
are modeled with relatively large disks. HD 100546 is modeled with a disk that flares more 
than that of HD 150193 and HD 179218. We find it interesting that the disk of HD 179218 
fiares less than the disk of HD 100546, even though the two are both Group I sources. At 
least on the quantitative level, the two disks differ. 

Although our model fits show no evidence for grain growth (all three SEDs are modeled 
with the same grain size distribution) there is a difference in the radial distribution of 
crystalline silicates. A possible evolutionary scenario is one in which when circumstellar 
disks form around young stars, they are thought to be primarly composed of amorphous 
ISM grains (Wooden, Woodward, & Harker 2004). As the disk evolves, crystalline silicates 
are condensed (Grossman 1972) or annealed (Rietmeijer et al. 2002) in the inner radial 
regions of the disk, either through heating (Brocklee-Morvan et al. 2002) or through shocks 
in the disk (Harker & Desch 2002). The crystals are then transported to the outer regions 
of the disk (Cuzzi, Dobrovolskis, & Champney 1993). Therefore, from this scenario, we can 
conclude that the disk around HD150193 is the least evolved since it does not contain any 
crystals. This is followed by HD179218 which has a large crystalline-to-amorphous in the 
inner radial region compared to the outer radial region. Finally, this is followed by HD 100546 
which has a slightly lower crystalline-to-amorphous silicate ratio at inner disk radii and a 
larger ratio at large disk radii compared to HD 1792 18. Such a scenario is supported by 
the findings of van Boekel et al. (2003) who find that the shape and strength (band over 
continuum) of the 10 /im silicate feature are correlated. Strong features (from small grains) 
show mostly amorphous type silicate grains, while weaker, fiatter features show evidence of 
crystaUine sihcates (grain processing). 

It should be noted that there may be evidence of grain growth in both HD 179218 and 
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HD 100546 as evidenced by the flatness of their 10 /xm sihcate features. We have not been 
able to adequately model "by eye" the shape of the feature with the current model. Attempts 
to use a more shallow grain size distribution and/or truncate the lower end of the grain sizes 
did not result in better fits (either by eye, or statistically) to the 10 fim feature. A more 
diverse mineralogy combined with a radially varying grain size distribution may improve the 
fits to the 10 /xm spectral feature. However, such models are beyond the scope of this paper. 



5.1. Comparison with Other Modeling Efforts 

Other authors have modeled the three objects presented here. As stated earlier, BdKDW 
modeled the SED of HD 100546 using a optically thin spherically symmetric dust distribution 
(i.e., no optically thick disk component). BdKDW find that the the amount of crystals used 
to fit the SED increased with radial distance from the star. This contrasts with our model 
which shows that there needs to be a higher crystalline-to-amorphous silicate ratio in closer 
to the star, in agreement with radial mixing models. BdKDW also finds that HD 100546 
has a large, flaring disk, consistent with our model results. 

Dominik et al. (2003), using another variation of the CG97 model (DuUemond, Do- 
minik, & Natta 2001) (hereafter DDN), used a simple mineralogy of grains 0.1 /im in size. 
They modeled HD 150193 with a disk smaller than that of the other two and found that 
the emission feature in HD 100546 is best flt with larger sized grains than the other two 
objects. While we can not conclusively make a quantitative statement about the grain size 
in HD 100546, our results are qualitatively similar with Dominik et al. (2003) for HD 150193. 
However, there are some important differences between our model results and those of Do- 
minik et al. (2003). They find that their fits for these three objects are best when the surface 
density increases with distance from the star (instead of decreasing as we have assumed here). 
They interpret these results as evidence for a "gap" in their disks, possibly produced by a 
planet or other large body. We do not find such evidence for an increasing surface density 
in our model calculations. Furthermore, based on DDN, it is extremely difficult to disentan- 
gle the degeneracy of the DDN model when simultaneously varying the surface density, the 
slope of the surface density with stellar distance, and the radial size of the disk. Therefore, 
while both of our results are intriguing, it underlines the need (as pointed out in Dominik) 
for high spatial resolution imaging of these objects in the mid-IR to resolve these modeling 
differences. 
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6. CONCLUSIONS 

In this paper wc constructed mid-IR SEDs of three Herbig Ae/Be stars using data from 
the HIFOGS, ISO SWS, IRAS, and MLOF photometry. We modeled the SEDs using the 
simple two-layer radiative transfer disk model of CG97 and COl with the added feature of 
calculating the emission from Mg-pure crystalline olivine grains. Our findings are highlighted 
as follows: 

1) HD 150193 has thermal emission arising from a relatively small (~ 5 AU in radius) 
flared {H/h = 3) disk. 

2) HD 179218 has thermal emission arising from a large (~ 150 AU in radius) less flared 
{H/h = 1) disk. 

3) HD 100546 has thermal emission arising from a large (~ 150 AU in radius) more 
flared {H/h = 4) disk. At 5 AU, the disk around HD 100546 flares 80% more than the disk 
around HD 179218 and 25% more than HD 150193. 

4) All three objects were modeled using a power law grain size distribution with a slope 
of a~^'^ for both the optically thick inner region and the optically thin surface layer. 

5) HD 150193 shows no evidence of emission from crystaUine sihcates while HD 179218 
and HD 100546 both show emission from Mg-rich crystalline olivine grains. 

6) The SED of HD 179218 was best flt using a crystaUine-to-amorphous ratio 84% larger 
in the inner radial regions of the disk (< 5 AU) compared to the outer regions. 

7) The SED of HD 100546 was best fit using a crystalline-to-amorphous ratio 37% larger 
in the inner radial regions of the disk (< 5 AU) compared to the outer regions. 

In circumstellar disks, grain growth and crystallization may require both sufficient disk 
mass (i.e., size) and scale height. Scale height may be a signature of turbulent movement 
in the disk which could aid grain growth and transport (Cuzzi, Dobrovolskis, & Champney 
1993). More disks need to be fitted with similar model structures to increase the sample to 
make more sense of the links between grain properties and disk structure. 

DEH and CEW acknowledge partial support for this work from from NSF Grant AST- 
0205814. In addition, DEH and DW acknowledge support from the NASA Ames Research 
Center Director's Discretionary Fund. The authors thank the day-crew and staff of NOAO 
CTIO for their assistance at the Blanco 4-m (in particular Ron Probst) and James E. Lyke 
for his assistance with the UM MLOF observations. Finally, the authors wish to thank an 
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Table 1. HIFOGS OBSERVATIONAL SUMMARY 





UT 


UT Time 


Integ. 


Airmass 


Flux 


Airmass of 


Scaling 


Object 


Date 


(hr:min) 


(min) 


Range 


Standard 


Standard 


Factor 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


HD 150193 


2000 May 13^ 


05:30 - 06:10 


20 


1.00 - 1.01 


7 Cru 


1.46 




HD 100546 


2000 May 13^ 


04:15 - 05:15 


32 


1.51 - 1.67 


7 Cru 


1.46 




HD 179218 


1998 June 23*^ 


07:03 - 07:23 


8 


1.15 - 1.11 


a Her 


1.14 


1.0 


HD 179218 


1998 June 24*^ 


07:32 - 08:10 


16 


1.11 - 1.12 


a. Boo 


1.63 


0.93 


HD 179218 


1998 June 25*' 


08:56 - 09:16 


8 


1.13 - 1.16 


/3Peg 


1.17 


1.0 


HD 179218 


1998 June 27^ 


07:52 - 08:00 


4 


1.11 - 1.11 


a Her 


1.18 


0.92 


HD 179218 


1998 June 30" 


05:24 06:00 


16 


1.29 1.19 


n Her 


1.13 


1.21 



^ Wavelength range: 7.4814 - 
Wavelength range: 7.5280 - 
Wavelength range: 7.8110 — 



13.222 //m 
13.122 /xm 
13.356 /xm 
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Table 2. IR PHOTOMETRY OF HD 179218 (2001 April 30.42 UT) 



Filter Observed Flux at Zero 

Wavelength Magnitude Magnitude 



(/xm) 




(W 


cm 


-2 




1.25 [J] 


7.19 ±0.10 


2, 


.90 


X 


10-" 


i.oo [Jri\ 


0. / 1 ± 0.09 


i. 




X 


1 n-13 
iU 


2.34 [K] 


5.88 ±0.03 


3 


.39 


X 


10-14 


3.65 [L] 


4.73 ±0.03 


6, 


.43 


X 


10-15 


4.90 [M] 


3.99 ±0.11 


1, 


.99 


X 


10-1^ 


10.00 [N] 


0.88 ± 0.03 


1 


.23 


X 


10-16 


7.91 


1.71 ± 0.06 


3 


.08 


X 


10-16 


8.81 


1.24 ± 0.05 


2, 


.02 


X 


10-16 


9.80 


0.81 ± 0.04 


1, 


.33 


X 


10-16 


10.27 


0.65 ±0.06 


1, 


.11 


X 


10-16 


11.70 


0.47 ±0.06 


6, 


.63 


X 


10-1^ 


12.49 


0.60 ±0.08 


5, 


.13 


X 


10-17 


18.30 


-1.00 ±0.09 


1, 


.14 


X 


10-17 


23.00 


-1.41 ±0.12 


4, 


.58 


X 


10-18 



Table 3. DATASET SCALING FACTORS TO HIFOGS 



Object 


ISO sws 


IRAS 


MLOF 


(1) 


(2) 


(3) 


(4) 


HD 150193 


0.90 


0.62 




HD 100546 


0.23 


0.22 




HD 179218 


1.08 


0.78 


1.0 
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Table 4. MODEL INPUT PARAMETERS 



Observed, 



Parameter 


Definition 


Fixed or Varied 


(1) 


(2) 


(3) 


MJMq) 


stellar mass 


observed ^ 




stellar radius 


observed ^ 




stellar temperature 


observed ^ 


= Ea3/2 (g cm-2) 


surface density at 1 AU 


fixed at 1000 


To (AU) 


outer disk radius 


varied 


H/h 


visible photospheric height/gas scale height 


varied 


Q 


slope of grain size distribution interior/surface 


varied 


fraax,s (A''-'^) 


maximum grain radius in surface 


fixed at 1 




maximum grain radius in interior 


fixed at 1000 


Tsub,Fe (K) 


iron sublimation temperature 


fixed at 2000 


Tsub,Ol (K) 


olivine sublimation temperature 


fixed at 1500 


T sub, ice (K) 


H2O ice sublimation temperature 


fixed at 150 



^ van den Ancker et al. (1998) 



Table 5. Minerals 



Dust lype 


Reference 


(1) 


(2) 


Amorphous Olivine 


Dorschner et al. (1995) 


Metallic Iron 


Pollack et al. (1994) 


Water Ice 


Warren (1984) 


Mg-pure Crystalline Ohvine 


Jager et al. (1998) 
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Table 6. Stellar Properties and Fitted Parameters 



Object 








To 


H/h 


Q 


f. a 
Ji,cr 


f 

J o,cr 


Mdisk 


H{ao)/ao 




(Me) 


(Re) 


(K) 


(AU) 










(Me) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


HD 150193 


2.3 


2.1 


9333 


5 


3 


3.5 






0.003 


0.11 


HD 100546 


2.5 


1.7 


10471 


150 


4 


2.5 


0.43 


0.27 


0.024 


0.42 


HD 179218 


4.3 


5.2 


10700 


150 


1 


3.5 


0.49 


0.08 


0.024 


0.09 



^ crystalline-to-amorphous silicate ratio for inner disk (< 5 AU). 

^ crystalline-to-amorphous silicate ratio for outer disk (5 - 150 AU). 

Disk mass (gas and dust) and maximum aspect ratio of the fitted disk are derived, not fitted 
or input, parameters. 
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Fig. 1.— HIFOGS spectra of HD 150193, HD 100546, and HD 179218. 
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Fig. 2.— Assembled SEDs for a) HD 150193, b) HD 100546, and c) HD 179218. Scaled to 
the HIFOGS spectrum {open circles) for each SED is the ISO SWS spectrum {black line), 
and the IRAS photometric points {gray squares) corresponding to each object. The 100 /xm 
point for HD 150193 (panel [a]) is an upper limit and not used in the model fitting. 
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CRYSTALLINE + 
AMORPHOUS 




Fig. 3. — Cross-sectional schematic of our disk model, adapted from Fig. 1 in COl. Identified 
in the sketch is the relative location of each mineral species. The vertical lines indicate the 
radial zones in the disk in which a particular mineral dominates the emission, and the dashed 
lines indicate the separation between the disk interior and the disk surface layer. The disk 
interior and disk surface layer {hash lines) are indicated. 
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Fig. 4. — COl model computed for three Herbig Ae stars: HD 150193, HD 100546 and 
HD179218. The disk interior (blue line), disk surface (red line), and stellar blackbody 
(orange line) are co-added to produce the model SED (black line). The model SED is 
compared to the assembled data sets including: ISO SWS spectra (indigo line), HIFOGS 
spectra (green circles) and IRAS photometry points (turquoise squares). HD 179218 also 
has MLOF photometry points (brown circles). 
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Fig. 5. — Zoom in on the COl model fits of tlie crystalline olivine features for the three 
Herbig stars: HD 150193, HD 100546 and HD179218. The colors are the same as in Fig. 4. 



